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Abstract 

A three-dimensional full Navier-Stokes (FNS) 
analysis was performed on a mixer/ejector nozzle 
designed to reduce the jet noise created at takeoff by a 
future supersonic transport. The PARC3D compu- 
tational fluid dynamics (CFD) code was used to study 
the flow field of the nozzle. The grid that was used in 
the analysis consisted of approximately 900 000 node 
points contained in eight grid blocks. Two nozzle 
configurations were studied: a constant area mixing 

section and a diverging mixing section. Data are 
presented for predictions of pressure, velocity, and total 
temperature distributions and for evaluations of internal 
performance and mixing effectiveness. The analysis 
provided good insight into the behavior of the flow. 

Nomenclature 

A area 

C D drag coefficient 

C fg gross thrust coefficient 

C p pressure coefficient 

L length of mixing section 

NPR nozzle pressure ratio 

P pressure 

T temperature 

V velocity 

W weight flow 

X distance from chute exit 

P mixing parameter based on velocity 

u mixing parameter based on total temperature 

Subscripts 

amb ambient 


bt boattail 

e nozzle exit 

mix mixing section 
o total 

p primary 

ref reference 

s secondary 

Introduction 

In the past decade, there has been an increasing 
interest in mixer/ejector nozzles because of the their 
noise suppression applications. These nozzles have the 
potential to lower jet noise without a significant thrust 
loss. 

The basic concept of forced mixer nozzles 1 ’ 6 was 
developed in early studies and showed the feasibility of 
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the nozzles. More recent work ’ has refined the designs 
and developed a practical system that could be 
employed on a future supersonic transport. Such an 
aircraft would have to meet Federal Aviation Regulation 
(FAR) 36 Stage III noise regulations. 

The mixer/ejector nozzles are designed to entrain 
large amounts of secondary flow through an array of 
lobed chutes that are deployed into the primary stream. 
These chutes alternately misalign the primary and 
secondary flow with respect to the axial direction 
(Fig. 1). The flow misalignment causes streamwise 
vortices to develop at each primary-secondary flow 
interface. These vortices rapidly mix the two streams 
together before exiting the nozzle. The rapid internal 
mixing of the two streams lowers the total jet velocity 
at the nozzle exit and results in lower noise; however, 
high thrust is maintained because of the large amount of 
flow augmentation. 

This paper examines a mixer/ejector nozzle at take- 
off conditions using a full Navier-Stokes (FNS) computa- 
tional fluid dynamics (CFD) code. The objectives of 
this study were to gain better insight into the nozzle 


flow field; to provide data for improvement of this 
design and of future nozzle designs; and to evaluate the 
ability of the CFD code to accurately predict 
mixer/ejector nozzle flow fields. 

Nozzle Description 

The General Electric (GE) Aircraft Engine Com- 
pany, under NASA Research Announcement (NRA) 
contract NAS-25415, has designed a two-dimensional 
(rectangular) mixer/ejector for noise suppression that 
could be used in conjunction with a mixed-flow turbofan 
engine. The nozzle consists of an array of 20 lobed 
chutes that deploy into the primary stream during 
takeoff and retract out of the flow path after takeoff. 
To avoid shocks in the flow field, the primary flow path, 
a convergent-divergent design, was designed to match 
the static pressure of the secondary stream when the two 
flows merge. The nozzle is intended to entrain approxi- 
mately 60 percent secondary flow. 

The nozzle studied in this paper (Fig. 2) is a 
simplified scale model of the aforementioned nozzle, and 
was designed to investigate the aerodynamic and mixing 
characteristics of the full-scale design. The general 
features remain the same, but the number of chutes was 
reduced from 20 to 10. The actuators, seals, hinges, and 
other components were removed to uncomplicate the 
model. GE’s Aerodynamics Research Lab (ARL) con- 
ducted an experimental investigation of this nozzle. 
Two nozzle configurations were analyzed in this compu- 
tational study. The baseline configuration has a con- 
stant area mixing section (A e /A mix = 1.0). The second 
configuration has a diverging mixing section 
(A e /A mix = 1.2). Figure 3 is a schematic of the nozzle 
showing terms used in the following sections. 

Analysis 

The PARC3D Code 

The PARC3D CFD code 9 was used in the analysis. 
The PARC code is a full Navier-Stokes multipurpose 
flow solver that was developed at the U.S. Air Force’s 
Arnold Engineering Development Center (AEDC). It 
is a central differencing code that utilizes a Beam- 
Warming approximate factorization algorithm. 10 One 
of the code’s distinguishing features is its ability to 
specify any portion of any grid surface as a boundary 
condition, which allows flexibility when generating a 
computational grid. Also, the code can solve grids that 
are made up of multiple grid blocks. Grid blocking 
allows one (l) to break up a complex domain into 
smaller sections, which separately use less computer 
memory than a single large grid and (2) to break up 
complex shapes into simpler components, which 
simplifies grid generation. Two grids that comprise a 


block interface do not have to be contiguous, which 
allows grid points across a block boundary to be 
reduced. A trilinear interpolation routine is used to 
transfer data between blocks. 

The code has several options to model turbulence. 
The default is an algebraic model based on the method 
developed by Thomas, 11 which can handle both wall- 
bounded and shear layer flows. Therefore, it can be 
used for a wide range of problems. A two equation k-e 
model, recently added to the code, is based on the work 
of Speziale and has been extended for use with com- 
pressible flow. The Baldwin-Lomax model, for wall- 
bounded flows is also available for use in conjunction 
with the Thomas model for flows with both walls and 
shear layers. The Thomas model was used for this 
calculation because of its simplicity and its speed. 

Grid Generation 

Because of the complexity of the nozzle geometry, 
grid generation was a very significant portion of the 
overall analysis. Figure 4 shows the computational 
domain used in this analysis. To reduce grid size and 
computational time, only one-half of a chute wavelength 
was modeled. Symmetry planes were specified on the 
primary and secondary flow centerlines; consequently, 
the effect of the nozzle sidewalls is not considered. In 
other words, the grid models an infinite row of chutes. 
The nozzle centerline was used as a plane of symmetry; 
therefore, it was necessary to model only the upper half 
of the nozzle. 

The grid used in the analysis consisted of 920 671 
grid points (Figs. 5 and 6). The grid was divided into 
eight smaller grid blocks to simplify grid generation. 
Each grid block is a component of the nozzle. For 
example, the chute, mixing section, and ejector inlet are 
all separate blocks of simpler shape than they would be 
if combined. Figure 6 shows an example of reduced grid 
points across a block boundary. The exterior flow 
upstream of the ejector inlet only requires resolution in 
the boundary layer. However, more resolution is needed 
near the leading edge of the ejector inlet. To avoid 
unnecessary high grid resolution upstream, a new coarser 
block was created in the upstream area, which saved 
approximately 67 000 points. 

The grid for each block was generated from detailed 
drawings on an IRIS workstation. The I3G/VIRGO 
interactive grid code 14 was used to define the geometry 
and create the six surface grids that comprised one grid 
block. The three-dimensional grid volume for the block 
was created by inputting the surfaces to the INGRID3D 
code. 15 The individual blocks were then combined in a 
post-processing step to form the complete flow-field grid. 
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Results and Discussion 


Results were obtained on both the Cray Y-MP and 
the Convex C220 computers at NASA Lewis Research 
Center. The Cray Y-MP calculated the solution at 
approximately 150 iterations per hr; the Convex, at 
11 iterations per hr. The Convex, although much slower 
in raw CPU speed, allowed the code to run contin- 
uously. The Cray Y-MP, on the other hand, used a 
queue system. Although the Cray was faster in terms of 
overall turnaround time, the difference between the two 
machines was not as significant as one would expect. 
These very low iteration speeds are due to the extremely 
large grid required to model the geometry. A case from 
an initial flow field to a converged solution requires 
approximately 40 000 iterations. Convergence was 
checked by monitoring the mass flow conservation 
through the nozzle. 

The data presented are for the following takeoff 
conditions; nozzle pressure ratio (NPR) of 4, primary 
stream total temperature of 850 °R, secondary total 
temperature of 530 °R, and a free-stream Mach number 
of 0.27. 

Nozzle Flow Field 

The Mach number distributions for the constant 
area mixing section (Figs. 7(a) and (b)) clearly show an 
overexpansion and resulting shock near the nozzle exit. 
The secondary centerline plot shows a region of higher 
speed flow beginning at the start of the shroud wall and 
extending downward toward the nozzle exit. Also, the 
stagnation point for the free-stream flow can be seen on 
the leading edge of the ejector inlet. 

The differences in total temperatures for the two 
streams distinguish them. Figure 8(a) shows that the 
temperature on the primary centerline does not decay 
until near the exit of the nozzle. However, the high 
temperatures appearing on the secondary centerline 
(Fig. 8(b)) indicates that some primary flow has been 
swept over into the secondary centerline plane by the 
vortex. This explains the region of high velocity in the 
mixing section (Fig. 7(b)). 

Figures 9 and 10 indicate a similar flow field for the 
A e /A m ix = 1-2 case. However, the flow separates from 
the shroud on the primary centerline at approximately 
X/L = 0.6. This large recirculating region does not ex- 
tend to the secondary centerline plane. 

Shroud Static Pressures 

Figure 11 shows the static pressure distributions on 
the shroud wall at both the primary and secondary flow 
centerlines. The areas of high pressure between 


X/L = 0.1 and X/L = 0.2 are due to a shock impinging 
on the shroud. Both nozzle configurations show the 
nozzle is highly overexpanded. For the constant area 
mixing section (A e /A mix = 1.0), the flow shocks at an 
axial location of 0.8 X/L to reach ambient pressure 
before exiting the nozzle. For the diverging mixing 
section (A e /A mix = 1.2), the shock caused by over- 
expansion occurs at different streamwise positions for 
the two spanwise locations plotted. The difference in 
shock location occurs because the recirculating region 
does not extend over the entire width of the nozzle. 

Primary and Secondary Flow Mixing 

Figure 12 is an illustration of the area for which 
quantities are plotted in Figs. 13 to 16. The most im- 
portant attribute of a mixer/ejector nozzle is its mixing 
effectiveness. The extent to which the two streams mix 
directly determines both noise characteristics and ejector 
performance. The enhanced mixing of the two streams 
in this nozzle is caused by the vortices created at the 
chute exits because of the misalignment of the primary 
and secondary flows (Fig. 13(a)). The ejector flow is 
drawn downward through the chutes. The primary 
stream flowing between the chutes is directed slightly 
upward. Downstream of the chute exit, the primary 
flow migrates upward and then begins to roll over and 
down into the secondary flow plane. The vortex forms 
near the top of the mixing section close to the shroud 
wall. As the flow moves through the nozzle, the center 
of the vortex translates downward and the size of the 
vortex grows. The motion of the higher momentum 
primary flow dominates and determines the motion of 
the vortex. DeJoode and Patankar show, in a hyper- 
mixing nozzle, that symmetric vortices are generated 
when the two flows have the same momentum. Fig- 
ure 13(b) also shows the formation of a small set of 
vortices at the shroud wall on the secondary flow center- 
line. These vortices persist throughout the nozzle. 

Total temperature contours (Fig. 14) are useful in 
visualizing the extent of mixing that has occurred. A 
mushroom-shaped plume of primary flow forms at the 
secondary flow centerline, downstream of the chute exit. 
This shape defines the vortex and the tip of the primary 
stream as it mixes with the secondary flow and is typical 
of mixer nozzle flow fields. At the nozzle exit, the 
primary flow extends down near the axis of symmetry. 
A considerable amount of mixing has taken place; 
however, significant amounts of unmixed primary flow 
that have migrated up near the shroud wall remain. 

For the diverging shroud, the flow behaves similarly 
to that of the constant area mixing section. However, 
the recirculating region at the nozzle exit forces the 
vortex away from the wall (Fig. 15). The total 
temperature contours (Fig. 16) show the recirculating 
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region at the exit, where the lower temperature ambient 
flow has been pulled into the region. The area extends 
over the entire primary flow path. 

2 

Bevilaqua proposed using the flatness of the 
velocity profile as a quantitative measure of the extent 
of mixing. This mixing parameter is defined as 



(V^A 


where the average velocity is given by 

<v> - (*> 

A 

For uniform velocity profile the parameter is /3 = 1.0. 
The unmixed flow at the chute exit plane has a parame- 
ter of /3 = 1.58. At the exit of the nozzle, the values for 
the A e /A mix = 1.0 and A e /A mix = 1.2 configurations 
are /3 = 1.19 and (3 = 1.30, respectively. Because the 
velocity field is dependent on more than just the extent 
of mixing, problems can arise when using this parameter 
to analyze the mixing of the two streams throughout the 
nozzle-mixing section. However, the total temperature 
field depends only on the mixing. In flows where the 
primary and secondary streams differ in total tempera- 
tures, this fact can be used to define a mixing 
parameter. 

Jl T o - T 0 JdA 

w s r? (3) 

J |T 0 - T |dA 
L J ave Jref 

The chute exit plane was used as the reference plane. A 
value of 1 for w corresponds to no mixing; a value of 
0, to complete mixing. Figure 17 shows the increase in 
mixing through the nozzle-mixing section. The mixing 
appears to occur at a constant rate. At the nozzle exit 
plane, the flows are 45- and 40-percent mixed for the 
constant area and diverging cases, respectively. Extra- 
polating from the data, complete mixing of the two 
streams would occur at an X/L of approximately 2.3. 
This indicates that if complete mixing is desired, very 
long and, consequently, very heavy nozzles would be 
required. 

Nozzle Performance 

Table 1 summarizes the major performance charac- 
teristics of the two configurations. The predicted ejector 
pumping corresponds well with the design value and is 
similar for both nozzles. However, the diverging shroud 


configuration was expected to pump larger amounts of 
flow. 1 Two possible causes for this behavior are (l) the 
recirculation region at the exit effectively reduces the 
nozzle exit area and (2) the secondary flow may have 
choked and the area increase downstream was not able 
to influence the pumping. 

Gross thrust performance is presented as a variation 
of the baseline configuration. The separated region in 
the diverging case caused a loss of thrust from the 

A e/ A mix = 10 case - 

Because of the large boattail angles at takeoff 
conditions, the pressure drag on the boattail may be a 
concern. Pressure coefficient on the boattail is plotted 
in Fig. 18. These curves are typical of flows over 
boattails. The boattail angle is larger for the 
A e /A mix = 1.0 case, which causes the lower surface 
pressures and, therefore, higher drag. 

Conclusions 

Mixer/ejector nozzles have the potential to lower jet 
noise without significant thrust loss. A full Navier- 
Stokes (FNS) analysis of a rectangular mixer/ejector 
nozzle was performed. The objective was to gain better 
insight into the complex flow field and to provide data 
for improvement of the design. 

The PARC3D code was used for the analysis. The 
grid that was used consisted of 920 671 grid points in 
8 grid blocks. The complex nozzle geometry required a 
large grid generation effort. Because of the large grid 
size, iteration speed was very slow. Two configurations 
of the nozzle, a constant area and a diverging mixing 
section, were studied at takeoff conditions. 

The flow field was dominated by a system of large 
streamwise vortices. These vortices were created at the 
exit of mixer/ejector chutes as a result of the misalign- 
ment of the primary and secondary flows. The vortices 
sweep the primary flow into the secondary stream, 
which increases the mixing between the two streams. 

The flow in both configurations overexpanded and 
shocked near the nozzle exit. For the diverging shroud 
case, the flow separated from the shroud and caused a 
large recirculating region which contributed to poorer 
thrust performance than occurred in the constant area 
mixing section. 
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Table 1. — Nozzle Performance 


Characteristic 

Mixing section configuration 

Constant area, 

A e /A mix = 1.0 

Diverging, 

A e /A mix = 1-2 

Mixing parameter, (5 

1.1985 

1.3111 

Mixing parameter, w 

0.5516 

0.6015 

Ejector pumping W 8 /W D 

0.6090 

0.6110 

Gross thrust coefficient, 

Baseline 

-6.74 percent 

Drag coefficient (boattail), ^D bt 

Baseline 

-99.56 percent 
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Figure 1 . — Typical mixer/ejector chute geometry. 



Chute Exit Plane — * 



Nozzle 

Centerline 
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Figure 3.— Nozzle Schematic. 



Figure 4. — Computational domain. 
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Figure 5. — Computational grid. 



(a) Primary centerline (Aq/A^ = 1 .0). 



(b) Secondary centerline (Aq/A^ = 1 .2). 
Figure 6. — Computational grid nozzle detail. 
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(a) Primary centerline. 



Figure 8. — Total temperature contours (Ag/A^x = 1 .0). 
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(a) Primary centerline. 



0.0 (b) Secondary centerline. 

Figure 9. — Mach number contours (Ag/A^ = 1 .2). 



(a) Primary centerline. 
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Axial position, X/L 

Figure 11 . — Shroud static pressures mixing section. 



Axis of symmetry 

Figure 1 2. — Area of mixing section. 
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(c) X/L = 1 .0. 

Figure 13.— Velocity vectors (A^A^ = 1 .0). 
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(a) X/L = 0.0. 



(c) X/L = 1.0. 

Figure 14. — Total temperature contours (A^/A^ = 1 .0). 


12 




I 


jWW'i m. 

WUllll///// 


m 1 1 


, Yilh .. 

\\\\ \ 1 1 1 / ////'i n*i 


1 1 1 1 1 hi 1 , 
ini iiip 
m 1 1 1 mi 

i i i i ini'"'" 1 1 ' ,//// mim i i i i i 1 1 mii NNN " 1 ' " //// imi 1 1 


\\\UI 1 1 1////' 1111 1 11111 1 Ml \\\ui | iu///' lln 1 1 

i i i i mi' 'ini 1 1 i i i 1 1 ini' 'in 1 1 1 

\\\\ 1111 / //// \\\\ 111 ////// 

'Mil I I I I I I I HIP 1 1 i||| | I | 


I ' I M l |||||.' ! "" l " , % im , , 

" v "" ' ' """'in i i i I I i i i iiiii"" ' ’ ' 1 mii 

• l I i l|‘ VSN ' ' ' ' ' 


I I I I I III' 

I I I » II lll v 


1 1 1* 


'll I I I I 
'''•III II I 


'll I I I I I I 
•It I I I I I I 


(a) X/L = 0.0. 


• " * * ' ' ’ r ^ , -.-V. >, -V \\\\\\ \ \ \ \ \ , # / / / ////////'/> 

— .^*N\\\\\\ Will// I ^_->>N\\\\\\ \ W I 11/ t !//////^_ 

/ ^vsW\\\\\ Will// Iff/S/s,^ ^.snWWW WII// tU////s,^ " K 1 

, wWWW WIN// //////✓.. % M J * . vsWWW Ml II// Iff///;,. . m M * 

j|,, vx\\miiiMi/////... M 1 J j , vwwu i 1 1 1 ;;/////, ^ | J J 

|j ,/ """""i""" j 1 1 /, — ' " i II M»/" ’ 

II 1 1 1 1 1 1 1 Ilium, ! | j/r .»uiiiilllllllllimm, M 


. ^^/tt/H H I I I I I 1 1 HHU\\ V ^ 


IW N 

I \ \ ' 


-.,,//! 1 I I II M \ Uw 


//mil II I ll\\ uuwn %< 


' y n i \ \ N 
* f / * 1 1 \ ' 


•///// / I I I I W\Wns 


''l l 

'"//I 


(b) X/L = 0.5. 

j j j |f rrrnm 


• i.uHMMUtU 1 I I I l/ll ////,,, . 

wvvWWI l 1 1 1 / ///////✓ 


nnnini 


»\\\\WWMI r I \ t I 1 1 ///,, , 
aWNUMII I II/// 


III'"' 

..^WUII"////^ 

*"M|||J / "' 

»nWU\M 1 / / / ///,,. ^ ^ 

III'- 

..sxxww w 1 1 tm,,,.. 

** M inii r ' 

««WW\\\I 1 II II//,, 

II/"*' 

>xww\\n 1 1 unn/,. 

"M \ 1 1 1 I I » / - * 
M|,i- 

.xwwwu l uni//,,.. % 

••ll 

III*" 

»wi\\umiirrf/i//. 

•'Minin'" 

*\mill|.lllllll/ln,. ni 

llll‘ l 

..//imiimimiiw. 

''miiiih"- 


l!U" 

^s//nti\i iiimnu,. 

'■'/IIIIU''- 

1 1 \ iiii 1 1 1 1 ^ j 

IIU' V 

—✓/// / // nmnwNNs^. 

n 1 1 1 1\"- 

— V/ f f | 


(c) X/L = 1.0. 

Figure 1 5.— Velocity vectors (Ag/A^ = 1 .2). 
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(a) X/L = 0.0. 



(b) X/L = 0.5. 



(c) X/L = 1 .0. 

Figure 16. — Total temperature contours (A e /A mjx = 1 .2). 


14 


Pressure coefficient, C p Mixing parameter, 



Figure 17. — Mixing parameter in mixing section. 



Axial position, X/L 

Figure 1 8. — Boattail pressure coefficient. 
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